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It is generally believed that the energy released in solar flares is stored 
Mh| in nonpotential magnetic flelds. The energy buildup could be a response 

O ' of the coronal magnetic field to the changes in photospheric magnetic 

environment caused mostly by sunspot motions and emerging fluxes. 
52 I Since the measurement of magnetic flelds in the corona is not available, 

the field measurement at the photospheric level has been widely used for 
the study of magnetic nonpotentiality in flare-producing active regions. 
/\ ' In particular, the change of some nonpotentiality parameters during 

3 . flare activity is regarded as a very important clue for understanding 

flare mechanisms. 

The MSEC (Marshal Space Flight Center) group studied the mag- 
netic nonpotentiality associated with solar flares using the MSEC 
magnetograph (Hagyard et al., 1982; Hagyard et al., 1984). They ex- 
amined the vertical current density, source field and magnetic shear 
derived from vector magnetograms of active regions in parallel with 
flare observations (Hagyard et al, 1984; Gary et al, 1987; Hagyard et 
ai, 1990). On the basis of these studies, they suggested a few important 
points characterizing flare producing active regions. They are a large 
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2 YONG-JAE MOON ET AL. 

shear angle, a strong transverse magnetic field along the neutral line 
and a twist of long neutral lines (Gary et al., 1991). The BBSO (Big 
Bear Solar Observatory) and Huairou (Beijing Observatory) groups 
have also presented several interesting studies on the change of non- 
potentiality parameters associated with solar flares (e.g., Wang et al., 
1996; Wang 1997). Wang (1992) defined a transverse field weighted 
angular shear and made a first attempt to study the change of the 
weighted mean shear angle just after a flare from vector magnetograph 
measurements. He found that the weighted mean shear angle jumped 
about 5 degrees, coinciding with a flare. He also showed that five more 
X-class and M-class flares had the same pattern of shear angle variation 
(see Table I in Wang 1997). It is to be noted that the shear increases 
in some active regions during flare activities coincide with emerging 
flux loops (e.g., Wang and Tang, 1993). On the other hand, the MSFC 
group reported that the shear may increase, decrease or remain the 
same after major flares (Hagyard, West, and Smith, 1993; Ambastha, 
Hagyard, and West, 1993). Also, Chen et al. (1994) observed that there 
were no detectable changes in magnetic shear after 18 M-class flares. 
Thus, the time variation of magnetic shear before and after a solar flare 
still remains controversial. 

For a quantitative understanding of nonpotentiality parameters such 
as angular shear, we have to keep in mind that there are several techni- 
cal problems calling for a careful treatment. First, the calibration from 
polarization signals to vector fields should be well established. Second, 
the 180° ambiguity of the observed transverse field has to be properly 
resolved. Most of nonpotentiality parameters such as vertical current 
density and magnetic shear critically depend on how well the 180° 
ambiguity problem is solved. Third, it is noted that Stokes line profiles 
could be changed due to the variation of thermodynamic parameters 
associated with flaring processes. 

So far, studies on the change of magnetic shear have been made by 
using mainly filter-based magnetographs thanks to their wide field of 
views and high time resolutions (Zirin, 1995; Wang et al., 1996; Wang, 
1997). However, filter-based magnetographs have some problems which 
originate from insufficient spectral information. Wang et al. (1996) 
well reviewed reliability and limitations of filter-based magnetographs. 
The calibration of filter-based magnetographs has often been made by 
employing the line slope method (Varsik, 1995) under the weak field ap- 
proximation (JefFeries and Mickey, 1991) or using nonlinear calibration 
curves based on theoretical atmospheric models (Hagyard et al, 1988; 
Sakurai et al, 1995; Kim, 1997). In these methods, the longitudinal 
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and transverse field components are independently given by 

Bl = kX , (1) 



Bj' = Kx 



I 

Q^ + U^ 



1/4 

(2) 



P 

in which Kl and Kt are calibration coefficients or curves which depend 
on the observing system or the models considered. Hagyard and Kineke 
(1995) developed an iterative method for Fe I 5250.22 line by consider- 
ing theoretical polarization signals with various inclination angles of the 
magnetic field. Recently, Moon, Park, and Yun (1999b) have devised 
an iterative calibration method for Fe I 6302.5 line following Hagyard 
and Kineke (1995). They applied this calibration method to the dipole 
model of Skumanich (1992) in determining the field configuration of a 
simple round sunspot. This study revealed that the conventional cal- 
ibration methods remarkably underestimate transverse field strengths 
in the inner penumbra. 

Most of previous studies on magnetic shear have adopted the po- 
tential method for resolving the 180° ambiguity (Wang et al, 1994). 
However, the method may possibly break down for fiaring active re- 
gions, which generally have a strong shear near the polarity inversion 
line. In filter-based magnetographs, magnetic field vectors are derived 
from polarization signals integrated over the filter passband and their 
calibration relations are determined without considering variation of 
thermodynamic parameters. That is to say, the derived magnetic fields 
may be affected by various physical processes involved in flares. 

In this study, we use a set of MSO magnetograms of AR 6919 
obtained by the Haleakala Stokes Polarimeter which provides simul- 
taneous full Stokes polarization profiles. Since the calibration is based 
on the nonlinear least square method (Skumanich and Lites, 1987), the 
magnetic field components derived are expected to be less affected by 
flare-related physical processes. Since it takes about an hour to scan 
a whole active region by the polarimeter, we can not detect such an 
abrupt change of magnetic shear as Wang (1992) found. In our study, 
emphasis is given to the evolution of magnetic nonpotentiality with 
the progress of an X-class flare which occurred in AR 6919. A similar 
study for AR 5747 is being prepared by Moon et al. (1999c). We expect 
that our data are able to yield a more reasonable estimate of various 
nonpotentiality parameters. In Section 2, an account is given of the 
magnetic nonpotentiality parameters that we have considered in this 
study. The observation and analysis are presented in Section 3 and 
the evolution of nonpotentiality parameters is discussed in Section 4. 
Finally, a brief summary and conclusion are delivered in Section 5. 



moon6919.tex; 1/02/2008; 11:31; p. 3 



4 YONG-JAE MOON ET AL. 

2. Magnetic Nonpotentiality Parameters 

2.1. Electric Current Density 

It is well accepted that electric currents play an important role in 
the process of energy buildup and relaxation of solar active regions. 
Since the observation of vector magnetic fields is available only in the 
photosphere, the vertical current at the photosphere is widely used in 
studies of solar active regions. According to Ampere's law, the vertical 
current density is given by 

in which fiQ = 0.012Gm/A is the magnetic permeability of free space. 

2.2. Magnetic Angular Shear 

Hagyard et al. (1984) defined the magnetic angular shear (or mag- 
netic shear) as the angular difference between the observed transverse 
field and the azimuth of the transverse component of the potential 
field which is computed employing the observed longitudinal field as a 
boundary condition. That is, the magnetic angular shear 9a is given by 

9a = 9o-ep, (4) 

in which 9o = &Tctan{By/Bx) is the azimuth of observed transverse 
field and 9p = arctan{Bpx / Bpy) is that of the corresponding potential 
field component. Noting that flares are associated with magnetic shear 
of strong transverse fields, Wang (1992) proposed a transverse weighted 
mean angular shear given by 

Oa = ^^, (5) 

in which Bt is the transverse field strength and the sum is taken over 
all the pixels in the considered region. In this work we use horizontal 
field strength in the heliographic coordinate instead of transverse field 
strength. 

2.3. Magnetic Shear Angle 

Lii et al. (1993) suggested a new nonpotentiality indicator, the angle 
between the observed magnetic field vector and the corresponding po- 
tential magnetic field vector. By definition, the shear angle 9s can be 
expressed as 

9s = arccos ( , ° n^^, I (6) 
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where Bq and Bp are the observed and potential magnetic field vectors, 
respectively. In this equation, Bp^ is identical with Bqz as explained 
earlier. In our study, we proceed a step further and consider a field 
strength weighted mean shear angle defined by 

where |B| is the field strength and the sum is taken over all the pixels 
in the considered region. 

2.4. Magnetic Free Energy Density 

The density of magnetic free energy is given by 

(Bq - Bp)^ Bs^ 
^^ 8^^ ^ ^^^ 

where Bg is the nonpotential part of the magnetic field, which was 
defined as the source field by Hagyard, Low, and Tandberg-Hanssen 
(1981). It can also be expressed as (Wang et a/., 1996) 

^^^ (|Bo[-jBp|)^ ^Mg^,in^(,^/2), (9) 

where |Bo| and |Bo| are magnitudes of the observed field and the 
computed potential field, respectively. The tensor virial theorem can 
be utilized to estimate the total magnetic free energy of a solar active 
region (e.g., Metcalf et al, 1995). However, as McClymont et al. (1997) 
pointed out (for details, see Appendix A of their paper), there are 
several controversial problems in estimating the magnetic free energy 
of a real active region with this method. In this study, we examine 
an observable quantity, the sum of magnetic free energy density over 
a field of view. This quantity is expected to indicate the degree of 
nonpotentiality at least near the photosphere. 



3. Observation and Analysis 

For the present work, we have selected a set of MSO magnetogram 
data of AR 6919 taken on Nov. 15, 1991. The magnetogram data were 
obtained by the Haleakala Stokes polarimeter (Mickey, 1985) which 
provides simultaneous Stokes I,Q,U,V profiles of the Fe I 6301.5, 6302.5 
A doublet. The observations were made by a rectangular raster scan 
with a pixel spacing of 2.8" (high resolution scan) or 5.6" (low resolution 
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Table I. Basic observational parameters of AR 6919. 



Data 


Date 


Time 


Scan Resolution 


Data Points 


AR 6919 a) 


Nov. 15, 1991 


17:50-18:55 


5.656" 


35x30 


AR 6919 b) 


Nov. 15, 1991 


21:05-22:45 


2.828" 


60x45 


AR 6919 c) 


Nov. 15, 1991 


23:46-25:25 


2.828" 


60x45 



scan) and a dispersion of 25 niA/pixel. Most of the analyzing proce- 
dure is well described in Canfield et al. (1993). To derive the magnetic 
fields from Stokes profiles, we have used a nonlinear least square fitting 
method (Skumanich and Lites, 1987) for fields stronger than 100 G 
and an integral method (Ronan, Mickey, and Orral 1987) for weaker 
fields. In that fitting, the Faraday rotation effect, which is one of the 
error sources for strong fields, is properly taken into account. The noise 
level in the original magnetogram is about 70 G for transverse fields 
and 10 G for longitudinal fields. The basic observational parameters of 
magnetograms used in this study are presented in Table I. 

To resolve the 180° ambiguity and to transform the image plane (lon- 
gitudinal and transverse) magnetograms to heliographic (vertical and 
horizontal) ones, we have adopted a multi-step method by Canfield et 
al. (1993) (for details, see the Appendix of their paper). In Steps 3 and 4 
of their method, they have chosen the orientation of the transverse field 
that minimizes the angle between neighboring field vectors and the field 
divergence |V-B|. Moon et al. (1999a) have already discussed in detail 
the problem of 180° ambiguity resolution for present magnetograms. 



4. Evolution of Magnetic Nonpotentiality 

Moon et al. (1999a) already explained in detail three magnetograms 
under consideration and their characteristics. Here we present the sec- 
ond image plane vector magnetogram of AR 6919 superposed on white 
light images in Figure 1. According to the Solar Geophysical Data, a 3B 
Ha fiare and an XI. 5 fiare occurred around 22:35 UT on November 15, 
1991 with the heliographic coordinate of S13 and W19. Considering 
the timing of the flaring events and the weak projection effect, our 
selected magnetograms are useful for studying the change of magnetic 
field structures associated with the fiare. This active region was also 
studied in terms of coordinated Mees/Yohkoh observations (Canfield 
et al., 1992; Wiilser et al., 1994), X-ray imaging observations (Sakao 
et al., 1992), white light flare observations (Hudson et al., 1992), and 
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Figure lA W?Pifiiage ptanS^lvfe© vecti^i^'lni.gfietogram of AR 6919 superposed on 
white light images taken on November 15, 1991. In the figure, the solid hnes stand 
for positive polarities and the dotted hnes for negative polarities. The contour levels 
correspond to 100, 200, 400, 800 and 1600 G, respectively. The length of arrows 
represents the magnitude of transverse field component. 



its preflare phenomena (Canfield and Reardon, 1998). From the three 
vector magnetograms in heUographic coordinates, we have derived a 
set of nonpotential parameters described in Section 2. In this study, we 
pay attention to their time variation associated the X-class flare. 

To examine the change of magnetic fluxes before and after the flare 
occurrence, we have plotted the time variation of magnetic fluxes of 
positive and negative polarities in Figure 2, in which the asterisked 
curves show the fluxes in a (5 sunspot region. As seen in the flgure, 
magnetic fluxes of the 5 sunspot region for both polarities decreased 
with time. It is also noted that there were several emerging fluxes 
(PI, P2, Nl, N2 in Figure 1) outside the 5 sunspot region, which were 
accompanied by eruptions of Ha arch fllaments (Canfleld and Reardon, 
1998). The new emerging positive fluxes such as PI and P2 account for 
total flux increase for positive polarities. These facts imply that the 
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Figure 2. Time variation of magnetic fluxes before and after an X-cIass flare in AR 
6919. The curves marked with an asterisk (*) stand for fluxes in a 5 sunspot region. 

Table II. Weighted mean magnetic shear angle, angular shear and sum of magnetic free energy den- 
sity for AR 6919 for three different times. All parameters are calculated for pixels with B^ > lOOG. 
The asterisked data are obtained adopting the potential field method for 180° ambiguity resolution. 





Data 
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22.3 


55.9 


33.5 


92.6 


3.6E3 


1.6E4 




1.1E23 


2.6E22 


c) 


23.0 
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X-class flare should be associated with both cancellation of the J-type 
magnetic flelds and new emerging fluxes. 

We show the vertical current density and inversion lines inferred 
from the longitudinal flelds in Figure 3. As seen in Figures 3-b) and 
3-c), the strongest vertical current density kernel is located near the 
inversion line of the 6 sunspot. As seen in the flgures, the vertical 
current density much increased just before the onset of the X-class 
flare. 
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Figure 3. Vertical current densities of AR 6919 superposed on white light images. 
The contour levels correspond to 3, 6, 9, 12 and 15 mA/m^, respectively. Thick solid 
lines represent the inversion lines inferred from the longitudinal magnetic field data. 
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a) Shear angle of AR6919 (91/11/15 17:50) 
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Position (arcsec) 

c) Shear angle of AR6919 (91/11/15 23:46) 




100 

Position (arcsec) 



Figure 4- Contours of angular shear multiplied by horizontal field strength obtained 
from three vector magnetograms of AR 6919 on Nov. 15, 1991 drawn in thick solid 
lines are superposed on the vertical magnetograms. (thick solid lines) The contour 
levels are 30000, 40000, 50000 and 60000 G deg, respectively. In all the figures, the 
thin solid lines stand for positive polarities and the thin dotted lines for negative 
polarities. The contour levels in the magnetograms correspond to 100, 200, 400, 800 
and 1600 G, respectively. 
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») Angular shear of AR6919 (91/11/15 17:50) 
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c) Angular shear of AR6919 (91/11/15 23:46) 
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Figure 5. Contours of shear angle multiplied by field strength drawn in thick solid 
lines are superposed on the vertical magnetograms. The contour levels are 30000, 
40000, 50000 and 60000 G deg, respectively. The magnetograms are the same as in 
Figure 4. 



moon6919.tex; 1/02/2008; 11:31; p. 11 



12 



YONG-JAE MOON ET AL. 



^^ 
o 

en 
QJ 

'So 



100 



80 



60 



40 



20 - 







Local 0^ 



Local 



18 



20 ?? 24 

Time (UT) 



26 



Figure 6. Time variation of weighted shear angle and angular shear before and after 
an X-class flare in AR 6919. In the figure Local 9a and Local 9^ are weighted shear 
angles for the S sunspot region. 



Figures 4 and 5 show angular shear weighted with transverse field 
strength and shear angle weighted with total field strength derived from 
vector magnetograms and computed potential fields (Sakurai, 1992). 
As seen in the figures, high values of shear angles are located near the 
S spot region. The changes of two mean shear angles with time are 
given in Figure 6, which shows that shear angles for the 6 spot area 
(denoted by Local in the figure) increased before the X-class flare and 
then decreased after it, but those for the whole active region decreased 
before the flare. It is to be noted that a similar pattern is also found in 
the evolution of a measure of magnetic field discontinuity, MAD (see 
Figure 7 in Moon et el. 1999a). The derived mean shear angles are listed 
in Table II, in which for comparison we also tabulate the corresponding 
values obtained with the potential field method for 180° ambiguity 
resolution. It is interesting that the mean shear angles for the 6 spot 
area (denoted by local in Table II) obtained with the potential method 
continuously increased during the fiaring activity unlike those obtained 
with our method. This eloquently demonstrates that totally different 
conclusions can be drawn on the shear angle evolution throughout a 
flare, critically depending on the 180° ambiguity resolution method 
employed in data reduction. 
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») Free energy density of AR6919 (91/11/15 17:50) 
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c) Free energy density of AR6919 (91/11/15 23:46) 
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Figure 7. Contours of free energy density obtained from three vector magnetograms 
of AR 6919 drawn in thick solid lines are superposed on the vertical magnetograms. 
The contour levels are 1 x 10'*, 2 x 10*, 3 x 10* and 4 x 10* erg/cm^, respectively. 
The magnetograms are the same as in Figure 4. 
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Figure 8. Time variation of the planar sum of free energy density before and after 
an X-class flare of AR 6919. In the figure Local corresponds to a 5 sunspot region. 



The magnetic free energy densities are shown in Figure 7. The free 
energy density near the 5 spot region largely decreased before the flare. 
We also list mean energy densities and total energy densities in Table 
II and plot their variations with time in Figure 8. These data indicate 
that magnetic free energy, at least in the neighborhood of the 5 spot, 
was indeed released during the flaring process. 

It is notably interesting and not easy to understand that in Figures 
6 and 8, the curvatures of curves for the whole active region are of 
opposite sign to those for the local 5 spot region. This may create 
some skepticism about the conventional picture of the flare mechanism. 
In Figure 6, where shear angle variations are displayed, the curves 
(Local 0a and Local 9s) for the 5 spot region seem to be consistent 
with the conventional flare picture, according to which the magnetic 
shear increases to build up free energy until the flare onset and then 
decreases as free energy being released. The free energy curve for the 
same 5 spot region in Figure 8, however, shows that the energy decrease 
started even before the flare onset although a more drastic energy 
release followed the onset. This indicates that relaxation of the twisted 
magnetic field mildly starts even hours before the visible flare onset. 
The decrease of mean shear angle for the whole active region before 
the flare supports this argument. Why, then, does the local shear in- 
crease before the flare? This may be interpreted as development of a 
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current sheet of a macroscopic size in the possible reconnection region 
near the 6 spot. Magnetic fields tend to take a lowest possible energy 
state under given constraints. The lowest energy state is expected to 
have a smooth spatial variation of magnetic field. However, a smooth 
state may not exist if the field has undergone too much twist. Then, 
a current sheet must develop in some regions while the field takes a 
smooth configuration in other regions. The process of current sheet 
development may be either a pure ideal MHD process as described 
above or a process involving a lot of small scale magnetic reconnections. 
The latter process can be analogized with an avalanche of a sand pile 
initiated by a slide of a few sand grains. This scenario was proposed 
and investigated by Lu et al. (1991, 1993). The increase of the planar 
sum of free energy density over the whole active region shown in Figure 
8 still requires more investigation. However, this can be also explained 
within the conventional picture of solar flares. When a magnetic field 
is twisted, the whole field tends to expand to occupy more volume per 
flux. The magnetic energy is thus more concentrated near the surface 
in a potential field than in a twisted field. The increase of the free 
energy after the flare onset might imply that the coronal field shrinks 
down while only part of the total free energy is released by the flaring 
process. However, our speculations based on only one flare case study 
definitely need to be examined by investigation of more cases. 



5. Summary and Conclusion 

In this study, we have studied the magnetic nonpotentiality of AR 6919 
associated with an X-class flare which occurred on November 15, 1991, 
using MSO magnetograms. The magnetogram data were obtained be- 
fore and after the X-class flare using the Haleakala Stokes Polarimeter. 
A nonlinear least square method was adopted to derive the magnetic 
field components from the observed Stokes profiles and a multi-step 
ambiguity solution method was employed to resolve the 180° ambiguity. 
From the 180° ambiguity-resolved vector magnetograms, we have de- 
rived a set of physical quantities characterizing the field configuration. 
The results emerging from this study can be summarized as follows: 

1) There was flux decreases for both polarities in a (5 sunspot pair as 
well as flux increases outside it, which implies that the energy release 
of the X-class flare should be associated with flux cancellation and 
emergence. 

2) It was found that the vertical current near the 6 sunspot much 
increased before the flare. 
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3) We also found that magnetic shear near the 6 sunspot increased 
before the flare and then decreased after it, while magnetic shear in the 
whole active region decreased before the flare. 

4) The sum of magnetic energy density decreased before the flare, 
indicating that magnetic free energy was released by the flaring event. 

However, we also found different evolutionary tendencies of nonpo- 
tentiality parameters for the whole active region and for the local 5 spot 
region. These differences must be examined through further studies 
of many other flare-bearing active regions. If some of them could be 
confirmed, these would serve as important clues to understand flaring 
mechanisms. 
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